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Abstract 

A linear dispersive mechanism for error focusing in polychromatic so- 
lutions is identified. This local error pile-up corresponds to the existence 
of spurious caustics, which are allowed by the dispersive nature of the 
numerical error. ^From the mathematical point of view, spurious caustics 
are related to extrema of the numerical group velocity. Several popular 
schemes are analyzed and are shown to admit spurious caustics. It isalso 
observed that caustic-free schemes can be defined, like the Crank-Nicolson 
scheme. 



1 Introduction 

The analysis and control of numerical error in discretized propagation-type equa- 
tions is of major importance for both theoretical analysis and practical applica- 
tions. A huge amount of works has been devoted to the analysis of the numerical 
errors, its dynamics and its influence on the computed solution (the reader is 
referred to classical books, among which ^0 [21 El)- It appears that existing 
works are mostly devoted to linear, one-dimensional numerical models, such as 
the linear advection equation 

du du 

where c is a constant uniform advection velocity. A striking observation 
is that, despite the tremendous efforts devoted to the analysis of numerical 



schemes in this simple case, the full exact non-homogeneous error equation has 
been derived only very recently 0]. 

The two sources of numerical error are the dispersive and dissipative prop- 
erties of the numerical scheme, which are very often investigated thanks to the 
spectral von Neumann analysis. Following this approach, a monochromatic wave 
is used to measure the accuracy of the scheme. Such a tool is very powerful and 
provides the user with a deep insight into the discretization errors. But some 
results coming from practical numerical experiments still remain unexplained, 
despite the linear character of the discrete numerical model. As an example, 
let us note the sudden growth of the numerical error for long range propagation 
reported by Zingg [7] for a large set of numerical schemes, including optimized 
numerical schemes. 

A limiting feature of the usual modal analysis is that it is applied to monochro- 
matic reference solutions. Therefore, dispersive phenomena associated to poly- 
chromatic solutions are usually not taken into account. 

The present paper deals with the analysis of a linear dispersive mechanism 
which results in local error focusing, i.e. to a sudden local error growth for 
polychromatic solutions. This phenomena is reminiscent of the physical one 
referred to as the caustic phenomenon in linear dispersive physical models 
and will be referred to as the spurious caustic phenomenon hereafter. It will be 
shown that, for some specific values of the Courant-Friedrichs-Lewy (number 
(CFL), spurious caustics can exist for some popular finite-difference schemes. 

The paper is organized as follows. The numerical schemes retained for the 
present analysis arc briefly recalled in section Main elements of cautic theory 
of interest for the present analysis are briefly recalled in sectional Test schemes 
are analyzed in section 0] 

2 Test numerical schemes 

For the sake of simplicity, the analysis will be restricted to three-level three- 
point numerical methods. The extension of the present analysis to other schemes 
is straightforward. A general finitc-diffcrcnccd version of the linear advection 
equation is 

au^+pu^+ju^ + Su^+vu^+eu^ + Cu^ + vu^ + eu^l = (2) 
with 

m m = u(lh,mr) (3) 

where h and r are the mesh size and time step respectively. For the sake of 
simplicity, these two quantities are assumed to be uniform. The CFL number 
is defined as a — cr/h, while the non-dimensional wave number is defined as 
ip = kh where k is the wave number of the signal under consideration. 
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A numerical scheme is specified by selecting appropriate values of the coef- 
ficients a, (3, 7, 5, e, £, r\ and 9 in Eq. J5J. Values corresponding to numerical 
schemes retained for the present works are given in Tabled 



Table 1: Numerical scheme coefficient. 
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3 Caustics 

The solution of Eq. Q is taken under the form: 

u(x,t,k) = e * (4) 

where u> = £ w + i 77^ is the complex phase, and fc the real wavcnumber. For 
dispersive waves, it is recalled that the group velocity V g (k) is defined as 



V B (k) - ^ (5) 



dk 

A caustic is defined as a focusing of different rays in a single location. The 
equivalent condition is that the group velocity exhibits an extremum, i.e. there 
exists at least one wave number k c such that 

§fe) = (6) 

The corresponding physical interpretation is that wave packets with charac- 
teristic wave numbers close to k c will pile-up after a finite time and will remain 
superimposed for a long time, resulting in the existence a region of high energy 
followed by a region with very low fluctuation level. 

The linear continuous model Eq. (y| is not dispersive if the convection ve- 
locity c is uniform, and therefore the exact solution does not exhibits caustics 
since the group velocity does not depends on k. The discrete solution associated 
with a given numerical scheme will admit spurious caustics, and therefore spu- 
rious local energy pile-up and local sudden growth of the error, if the discrete 
dispersion relation is such that the condition (jSJ) is satisfied. For a uniform 
scale-dependent convection velocity, such spurious caustics can exist in poly- 
chromatic solutions only, since they are associated to the superposition of wave 
packets with different characteristic wave numbers. 
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The general dispersion relation associated with the discrete scheme is 



ae l * + ^ e lip + C] e 2tulT + {{(3 + e 11 ? 5) + e}e lv e lu:T +rie 2tlp + 6 = (7) 

which is a non linear quadratic equation in e luJ T that can easily been solved. 
The corresponding group velocity is given by: 

(8) 

The numerical solution will therefore admits spurious caustics if 

dk dip dk dip 

The corresponding values of ip and k will be respectively denoted ip c and k c . 
Spurious caustics are associated with characteristic lines given by 



where 



~ = Uc (10) 



U c = V g {p c ) (11) 



We now illustrate the caustic phenomenon considering the two following 
sinusoidal wave packets: 



ui = e -°t*-*l-°*)' Cos[fci (x-xl-ct)] , u 2 = e -<*(x~*4-ct) 2 Cos[fc 2 {x-xl-ct)} 

(12) 

where a > 0). The two wave packets are initially centered at x\ and Xq, 
respectively. 

If the solution obeys the linear advection law given by Eq. , the initial field 
is passively advected at speed c, while, if the advection speed is scale-dependent 
(as in numerical solutions), the two packets will travel at different speeds, leading 
to the rise of discrepancies with the constant-speed solution. Another dispersive 
error is the shape-deformation phenomenon: due to numerical errors, the exact 
shape of the wave packets will not be exactly preserved. This secondary effect 
will not be considered below, since it is not related to the existence of spurious 
caustics. 

The spurious caustic will appear if the two wave packets happen get su- 
perimposed. During the cross-over, the L m norm of the error (defined as the 
difference between the constant-speed solution and the dispersive one) will ex- 
hibit a maximum. The characteristic life time of the caustic depends directly 
on the difference between the advection speeds of the two wave packets and the 
wave packet widths. 
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Neglecting shape-deformation effects and assuming that the numerical scheme 
is non-dissipative, the numerical error E is given by: 



E = | e- a ^- x o- c * )2 Cos[/ci (x - xl - ct)} - e -«(*-*o-* V0 a Cos[fci(a: -x\- tVi)] 

+ e - a{x - x »- ct)2 Cos[k 2 (x - x 2 - ct)] - e -«(*-*S-*v a ) a Cos [fc 2 (x -xl-tV 2 )] 

(13) 

where Vi and V% are the advection velocity of the two wave packets, respec- 
tively. 

A simple analysis show that 



lim £«,(£(<)) =ioo(«i(* = 0)), maxL oo (£(t))=2L oo (u 1 (t = 0)) (14) 

i — >+oc i 

The time histories of the L\ and norms of E are displayed in Fig. ^ 
showing the occurance of the caustic and the sudden growth of the L m error 
norm. 
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Figure 1 : Time history of the numerical error for the two- wave packet problem 
(shape deformation and dissipative errors are neglected to emphasis the linear 
focusing phenomenon) . Top: L\ norm. Bottom: norm. Numerical parame- 
ters are a = 1, h = 0.01, V\ = 2.04, V 2 = 2.02, corresponding to the properties 
of the Lax scheme for a = 0.7 
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4 Analysis of test numerical schemes 

We analyzed the properties of the numerical scheme displayed in section [5] 
For the Leapfrog scheme, the dispersion relation is: 

l-e 2 ^ r + 2icrSin \p\e iur = 
from which it comes that 

£ w r = ArcSin [a <p] 
The group velocity can be expressed as 

c Cos [ip] 



y 9 = i 

yl-<r 2 Sin 2 [<p] 

leading to 

dV g _ ch3m[<p] ( ct 2 Cos 2 



- 1 



d * y/l-^Sm 2 MU 1 -^ Sin2 M) 
A trivial root is ip = mod tt, which corresponds to 

U c = ±c 

Zeros of Eq. (|19|l in the (u, (p) plane are displayed in Fig. [3 
Let us now consider the Lax scheme. The dispersion relation is: 

2e lv + {{a - 1) e 2 v - {a + 1)} e luJT = 

which yields: 



e^ T = -y/Cos 2 [<p] + a 2 Sin 2 [<p] 

t^Lj t = ArcTan [ a Tan [<p] ] 
The group velocity is equal to 

c 

9 ~ Cos 2 [</?]+ a 2 Sin 2 M 

yielding 

dV g 2r 3 (a 2 - 1) Sin[<p] Cos [<p] 



9 ^ a (Cos'^+^Sin 2 ^])' 



<p 
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Figure 2: Roots of Eq. I|19fl . which correspond to possible occurance of spurious 
caustics for the Leapfrog scheme. 

Spurious caustics (see Fig. [21) can arise for ip = mod it and ip — ^ mod ir. 
The corresponding group velocities are respectively: 



U c L c, U c z ^ 

G 

For the Lax-Wendroff scheme, the dispersion relation is: 

v +a(a- l)e 2iv - a {a + 1) + 2 (cr 2 - l)e lip e luJT =0 
which yields: 



1-4(1 -a 2 ) a 2 Sin 4 [|] 



c Sin [ip] 



I -2a 2 Sin 2 [f] 
The group velocity is given by 

{ (1 - 2a 2 Sin 2 [f ] ) Cos [cp] + a 2 Sin [cp] Sin 2 [f ] } 



V g = c- 



;i-2a 2 Sin 2 [§]) +a 2 Sin[(^] 2 



(26) 
(27) 

(28) 
(29) 

(30) 



7 



0.2 0.4 D . 



Figure 3: Roots of Ea. l|25|) . which correspond to possible occurance of spurious 
caustics for the Lax scheme. 

from which 



dV g c 2 t (-1 + a 2 ) (-2 + a 2 + 3 a 4 - 4 a 4 Cos [<p] + a 2 (-1 + a 2 ) Cos [2 <p]) Sin fo>] 



dip 



2 cr ( (1 - a 2 + a 2 Cos [p]f + a 2 Sin 2 [<p] 



(31) 

Roots of Eq. I|31|l arc shown in Fig. 

We finally consider the Crank- Nicolson scheme, whose dispersion relation is: 



a 2 ( Cos [<p] - 1) Cos [— ] + i cr Sin [— ] = 



which yields: 



„,„ T= ±(cr-2a 2 Sin 2 [f]) 
C7 ^ 2 cr 2 Sin 2 [f] 

(„,t = mod it 



(32) 



(33) 



(34) 

The group velocity being constant, no spurious caustics can arise with this 
scheme. 



5 Concluding remarks 



The existence of spurious numerical caustics in linear advection schemes has 
been proved. This linear dispersive phenomenon gives rise to a sudden growth of 
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Figure 4: Roots of H^ljl . which correspond to possible occurance of spurious 
caustics for the Lax-Wendroff scheme. 

the Loo norm of the error, which corresponds to a local focusing of the numerical 
error in both space and time. In the present analysis, spurious caustics have 
been shown to arise in polychromatic solutions. 

The energy of the caustic phenomenon depends on the number of spectral 
modes that will get superimposed at the same time. As a consequence, the 
spurious error pile-up will be more pronounced in simulations with very small 
wave-number increments. 

It has been shown that most popular existing schemes allow the existence of 
spurious caustics, while some schemes are caustic-free, like the Crank-Nicolson 
scheme. 
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